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ABSTRACT 

Cosmic structure formation is characterized by the complex interplay between gravity, turbulence, 
and magnetic fields. The processes by which gravitational energy is converted into turbulent and 
magnetic energies, however, remain poorly understood. Here, we show with high-resolution, adaptive- 
mesh simulations that MHD turbulence is efficiently driven by extracting energy from the gravita- 
tional potential during the collapse of a dense gas cloud. Compressible motions generated during 
the contraction are converted into solenoidal, turbulent motions, leading to a natural energy ratio of 
Esoi/Etot ~ 2/3. We find that the energy injection scale of gravity-driven turbulence is close to the 
local Jeans scale. If small seeds of the magnetic field are present, they are amplified exponentially fast 
via the small-scale dynamo process. The magnetic field grows most efficiently on the smallest scales, 
for which the stretching, twisting, and folding of field lines, and the turbulent vortices are sufficiently 
resolved. We find that this scale corresponds to about 30 grid cells in the simulations. Wo thus 
suggest a new minimum resolution criterion of 30 cells per Jeans length in (magneto)hydrodynamical 
simulations of self-gravitating gas, in order to resolve turbulence on the Jeans scale, and to capture 
minimum dynamo amplification of the magnetic field. Due to numerical diffusion, however, any exist- 
ing simulation today can at best provide lower limits on the physical growth rates. We conclude that 
a small, initial magnetic field can grow to dynamically important strength on time scales significantly 
shorter than the free-fall time of the cloud. 

Subject headings: hydrodynamics — magnetohydrodynamics — dynamo — ISM: clouds — ISM: kine- 
matics and dynamics — ISM: structure — methods: numerical — turbulence 



1. INTRODUCTION 

Astrophysical fluids on virtually all scales are charac- 
terized by highly complex turbulent motions. Besides be- 
ing turbulent, we know from observations that the gas is 
magnetized between galaxies and inside galaxies, as well 
as in individual star-forming molecular clouds, down to 
the proto-stellar accretion disks, which naturally accom- 
pany the birth of stars and planetary systems. However, 
not just the present-day universe is characterized by tur- 
bulence and magnetic fields. With the advent of large- 
scale computer simulations, it became clear that primor- 
dial gas becomes highly turbulent during the formation 
of the first galaxies and stars as well (e. g., Abel et al. 
2002l |O'Shea fc Nor man|[20WllWise fcXbel 2007, Pa^ 



et al. 2008 ; Grcif ct al. 2008 ) , which has far-reaching con 
sequences for all cosmic structure formation. This turbu 
fence provides an extremely efficient wa y to amplify pri 
mordial seeds of the magnetic field (e.g. , 
■Rulsrud ct al. 19971 [Lazar et al |P059{ 



B ierm ann 1950 

, ~t al. 2005 

Miniati fc: Bell||2011) by the turbulent dynamo process 



see the review by Brandenburg & Subramanian 2005 1 
already during the earliest phase s of structure formation 
in the universe ( Ryu et al.|2008 ). Here the amplification 
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of the magnetic field arises from sequences of a "stretch, 
twist, and fold" nature, until the magnetic field lines are 
so tightly packed that the magnetic energy density be- 
comes comparable to the kinetic energy density. Strong 
magnetic fields, close to equipartition may have been 
reached in normal disk galaxies already at redshif ts of 



z < 3 (Bernet et al. 2008 Arshakian et al. 2009 and 



references therein), emphasizing the significance of the 
dynamo action even in the early stages of galaxy evolu- 
tion and star formation. Thus, magnetic fields may not 
only have a significant impact on the gas dy namics and 



fragmentation in present-day star formation ( Hennebelle 
fc: Teyssi"er||2008 ), but potentially also in priniordial star 



formation. 

Despite it s ubiquity and importanc e for star formation 
( Mac Low fc Klesscn 2004 ; Elmcgrccn fc Scalop004; Mc] 



Kee fe Ostr ikcr 2007j), only very little is known about the 
origin of astrophysical t urbulence and mag netic fields 



Recently , Klessen fc He nnebelle (2010), Elmegreen fc 
Burkert (2010l), and Vazquez-Semadeni| ( )2010( ) suggested 



that a universal source tor driving astrophysical turbu 
lence everywhere from cosmic scales and galaxies, over 
galactic clouds, down to individual proto-stellar disks is 
the accretion from their environment, i.e., by gravity - 



driven flows, on all scales (see also. Field et al. 



Similar ideas existed earlier. In particular, Hoyle 



2008). 



showed that during the contraction of a slightly unsta- 
ble, nearly isotherm al gas cloud, internal ra ndom mo 



tions can be excited. Odgers fc Stewart ( 1958") suggested 
that even fully irrotational motions can give rise to tur- 
bulence due to the Reynolds stresses. Hoyle' s idea of 
gravity-driven tu rbulence was later refined by |Scalo fc| 
Pumphrey (1982), which they called 'turbulent virializa- 
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tion'. Fleck (1983) suggested that the injection of tur- 
bulence by gravitational contraction is important in the 
interstellar medium, and Begelman & Shlosman (20091 
conclude that angular-momentum transport during the 
turbulent collapse of a gaseous system may suppress frag- 
mentation. 

Here we show with numerical simulations that turbu- 
lence and magnetic field growth are indeed efficiently 
driven by the gravitational energy released during the 
collapse of a dense gas cloud. The connection between 
gravity-driven turbulence and magnetic field amplifica- 
tion has also been su ggested recently in a mode l by 



Schleicher et aL| ("2010[) and confirmed numerically in Sur 
In 



et al. 1 2bl0p . in this process, potential energy is con- 



verted into turbulent motions, which in turn amplify the 
magnetic energy via the turbulent dynamo. Thus, the 
driving of turbulence and magnetic field growth by gravi- 
tational infall may be the consequence of a self-regulating 
instability. At the bottom of this c ascade, close to the 
sonic scale (Federrath et al. 2010b), gas is expected to 



become subsonic as a consequence of a steep rise in the 
temperature when the gas becomes optically thick. Dur- 
ing this process, compressible modes will be converted 
into solenoidal turbulent motions, u ntil a natural energy 
ratio of Eso x/Etot ~ 2/3 is re ached, (Elmegreen & Scalo 



2004 Federrath et al. 2008b) 



In this paper, we first discuss the physics of magnetic 
field amplification by gravity-driven turbulence, and sub- 
sequently derive a new resolution criterion required to 
resolve these processes. We aim to address the follow- 
ing key questions: On what scales does the magnetic 
field grow during the collapse of a dense, magnetized 
gas cloud? What is the effective kinetic energy injection 
scale of gravity-driven turbulence? How are the com- 
pressible motions in a contracting system converted into 
turbulent random motions, and what is the asymptotic 
fraction of solenoidal motions generated during the con- 
traction? After presenting our methods in section [2] we 
address these physical questions in section [3] with the 
use of magnetohydrodynamical (MHD) simulations. We 
find that the magnetic field is most efficiently amplified 
on the smallest resolvable scales in the simulations, and 
grows exponentially fast due to the small-scale dynamo. 
The effective energy injection scale of gravity-driven tur- 
bulence is close to the local Jeans sc ale during the con- 
traction. Finally, we show in section 4.2 that about 2/3 
of the total kinetic energy released during the collapse is 
converted into solenoidal, turbulent motions, efficiently 
driving magnetic field amplification. 

In the second part of the paper (section |4]) , we dis- 
cuss the numerical resolution requirements for modeling 
turbulent, self-gravitating systems and for minimum dy- 
namo action to set in. To study the resolution depen- 
dence of our results, we use a sequence of simulations 
in which we resolve the Jeans length w ith 8, 16, 32, 64, 
and 128 grid cells (see, Sur et al^|2010 paper I). With a 
Fourier analysis of the magnetic energy, we confirm our 
earlier findings in paper I, showing that 30 grid cells per 
Jeans length is the threshold resolution for minimum dy- 
namo amplification of the magnetic field to occur. In 
addition, we show that a Jeans resolution of about 30 
grid cells is required to obtain converged values of the 
turbulent energy, in particular of the solenoidal (rota- 
tional) component, which drives dynamo amplification. 



In contrast, we find that the turbulent energy density, 
i.e., the turbulent pressure on the Jeans scale is under- 
estimated, if the Jeans length is resolved with 16 grid 
cells or les s during the colla pse. Apart from a few excep- 
Abel et al.|2002 ), the Jeans length is resolved 
To" 



tions (e.g. 



with less than 16 grid cells in almost all numerical sim- 
ulations today. We speculate th at this is because in the 
study by Truelove et al. (1997), it was found that only 
four grid cells per Jeans length are enough to avoid ar- 
tificial fragmentation. An equivalent resolution criterion 
for the Jeans mass in smoothed particle hydrodyn amic s 
(SPH) simulations was found by [Bate fc Burkert] fl997| ) . 
Also, the computational expenses increase strongly, if one 
aims to resolve the Jeans length with more than a few 
cells (or few particles in SPH). Thus, most existing hy- 
drodynamical and MHD simulations of self-gravitating 
media have typically used about ten grid cells per Jeans 
length or less. Some modification of this criterion was 



recently also suggested by Gawryszczak et al. (2010) to 
better resolve self-gravitating disks. 

Here we find that the turbulence and magnetic field 
structure are under-resolved, if the Jeans length is re- 
solved with 16 cells or less. Moreover, turbulent dynamo 
amplification of the magnetic field is completely absent in 
this case. To avoid this problem, we suggest a new reso- 
lution criterion for simulations of self-gravitating gaseous 
media: in order to resolve turbulence on the Jeans scale 
and to account for the turbulent pressure on that scale, 
as well as to capture minimum dynamo action in MHD 
simulations, 30 grid cells per Jeans length are required, 
which is discussed in detail in section |4l 

2. METHOD 

Gravitational collapse, turbulence and magnetic field 
evolution are -except for some idealized cases- diffi- 
cult to study analytically, because the system is highly 
non-linear and naturally three-dimensional. Thus, we 
study the processes leading to turbulence and magnetic 
field amplification through gravitational collapse in high- 
resolution MHD simulations. 

2.1. Initial conditions and setup of our MHD 
simulations of a collapsing dense core 

We present results of a numerical experime nt with a 



collap sing, magnetized, turbulent gas core (see, Sur et al 



2010[ paper I). We focus on the gravitational collapse and 



magnetic field amplification of a dense gas cloud, using a 
simplified setup, where we assume an almost isothermal 
equation of state (effective F — d\ogT/d\ogp+ 1 = 1.1, 
with the temperature T and density p) and neglect non- 
ideal MHD effects (discussed below). The numerical sim- 
ulations presented here were performed with the pub- 
licly availa ble adaptive-mesh refinement (AMR) code, 
FLASH2.5 (iFryxell et al.|[2000l ). We solve the equations 
of ideal MHd, including selt-gravity with a refinement 
criterion guaranteeing that the Jeans length. 



7rc£ 
Gp 



1/2 



(1) 



with the sound speed Cg, gravitational constant G, and 
density p, is always resolved with a user-defined num- 
ber of cells. We use the new HLL3R scheme for ideal 
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MHD ( Waagan et al.|20TT I , which employs a 3- wave ap 



proximate MiiD Kiemann solver ( Bouchut et al. 2007 



Waaga n"2009! i Bouchut et "al]|2010 ). This MH U scEeSe 



IS an ex tension ot the hydrodynamical version (Klingen 
berg et al.||2007[ ), developed for FLASH that preserves 



physical states (e.g., positivity of mass density and pres- 
sure) by construction, and is highly efficient and accurate 
in modeling astroph ysical MHD problem s involving tur- 
bulence and shocks ( Waagan et al.|[2"011 ). 

The ionization degree is expected to be sufficiently high 
in primordial clouds to ensure a strong coupling betw een 
ions and neutrals to maintain flux-freezing pvlaki fc Susa| 



2004 1. However, non-ideal MHD effects may eventually 
become important at very high densities, as suggest ed by 
simulations of contemporary star formation (e. g 



nebelle fc Teyssier 
ettects are not inc. 



20081 puffin Pudritz|[2009 | 



Hen 

'hese 



uded in the present calculations, but 
should be subject of future studies. 

The initial conditions for our numerical simulations 
were motivated from l a rge-scale cosmologica l models 



(e.g., Abel et al. 2002 Yoshida et al. 2008). How- 



ever, for simplicity we use a spherically symmetric initial 
density distrib ution, follow i ng a super-cri tical Bonnor 



Ebert profile (Ebert 1955 Bonnor 1956) with a core 
density of pbe = 4.68 x 10~^^ gem""* at a temperature 
of T = 300 K. We did not include any initial rotation 
of the core. The Bonnor-Ebert sphere has a radius of 
1.5 pc and a total mass of M = 1500 M©. The three- 
dimensional Cartesian computational domai n has a vol 



ume of (3.9 p c) . For a similar approach see Clark et al 



( |2008] |2011[ ) 

Although the existence of magnetic fields on a wide 
range of astrophysical scales is clear, their origin i s still 
poorly under stood. Some mechanisms include the Bier- 



mann 



et al. 



( 1950 ) battery, the Weibel instability (e.g., ^Lazarj 
2009p , the ejection of magnet i c fiel ds from ac- 



tive galactic nuclei (e.g., Xu et al. 2009[ and refer 



ences the rein), production by cosmic rays (Miniati & 



[BeU"2011), and amplification by galaxy cluster mergers. 
Recent lower limits of magnetic fields i n the intergalac- 
tic inedium suggest -Brm s ^ 10~^^ G (|Tavecchio et al 



2010[ [Dolag et al.||2011[ ). To keep our setup as simple 
and general as possible, we add a small initial seed of 
the magnetic field and study its self-consistent amplifi- 
cation during the collapse of a gravitationally unstable 
gas cloud in detail. In order to initialize seeds of tur- 
bulence and magnetic field, we impose a random initial 
velocity field with transonic velocity dispersion of am- 
plitude 1.1 km s""'^ (equal to the initial sound speed) and 
a weak random magnetic field with i?rms — 1 nG, with- 
out a mean field component. This initial field strength 
corresponds to a plasma beta (ratio of thermal to mag- 
netic pressure) of /3 « 10^". Both, the turbulent velocity 
and magnetic field were set up with the same power-law 
dependence in wavenumber space, P{k) cx fc~^, peaking 
on scales of about 0.8 pc, which roughly corresponds to 
the initial Jeans length of the core. The choice of the 
initial power spectrum is motivated by the finding that 
turbulence develops a k~^/^ and fc~^ v elocity spectrum 



for incom pressible Kolmogorov ( 1941 ) and for shock- 
dominated Burgersf l 1948^) turbulence, respectively. As- 
trophysical sysEems typically exhibit spectra of the tur- 
bulence in between these two extremes. Since both, Kol- 



mogorov and Burgers turbulent spectra are dominated 
by large-scale modes, i.e., small wavenumbers fc, our re- 
sults are not expected to depend strongly on the par- 
ticular choice of the initial spectrum. However, varying 
parameters of the initial conditions will be subject of a 
follow-up study. After an initial transient phase, the tur- 
bulence develops its own self-consistent spectrum, driven 
by the gravitational collapse. 

We consider an isolated system for solving the Poisson 
equation, and use Neumann boundary conditions (zero- 
gradient) for the MHD. Since the initial conditions for 
the magnetic field were generated in Fourier space to ob- 
tain a divergence-free field, the initial magnetic field is 
naturally periodic. This introduces an inconsistency at 
the six boundaries of the simulation domain and may 
cause some non-vanishing divergence there. However, 
even the fastest MHD waves (the fast magneto-sonic 
waves) have a box-crossing time of about 3 x 10^ yr, an 
order of magnitude longer than the initial free-fall time 
scale of the system. Thus, any potential perturbations 
produced at the boundaries can be ignored in the dense 
core. 



C onsistent with pr e vious works of Omukai et al. ( 2005 1 
and Glover & Savin (2009), we adopt an effective equa- 
tion oT'state^witirT = 1.1 for number densities in the 
range n — 10^-10^" cm^'^. However, this kind of setup is 
sufficiently general to allow us to investigate basic prop- 
erties of turbulence and magnetic field amplification in 
the context of collapsing primordial as well as present- 
day gas clouds. 

2.2. Five simulations with increasing Jeans resolution 

We note that the efficiency of the turbulent dynamo de- 
pends on the kinematic Reynolds number. Re = (intO'v/i^, 
with the driving scale £i^t, turbulent velocity dispersion 
(T„ on that length scale, and physical viscosity j/. Also 
note that previous studies found a critical value of the 
magnetic Reynolds number, Rm = tmtO' v/v, below which 
the sm all-scale dynamo is not excited (Brandenburg & 
Subramanian 2005). In numerical simulations, however, 
the viscosity, ;/ and the magnetic diffusivity, rj are typi- 
cally determined by the numerical cutoff scale, as a result 
of the discretization of the fluid variables on a compu- 
tational grid with finite resolution. The numerical val- 
ues of Re and Rm are thus related to the grid resolu- 
tion and to how well the turbulent motions are resolved 



(e.g.,|Haugen et al.[2004[ Balsara et al.||2004[ |Schekochi 
hm et al. 120041). Even if physical viscosity /dittusivit-\ 



hin et al.||2004p . Even if physical viscosity/ dittusivity 
IS added as an extra term in the equations, it must be 
guaranteed that it is at least as high as the numerical 
viscosity/diffusivity. Typical Reynolds numbers reached 
in numerical simulations of turbulence are of the order of 
a few hundred (depending on what scales are estimated 
to be affected by nu merical viscosity, whic h varies among 



different codes, see Kitsionas et al. 2009), which is usu- 



ally orders of magnitude below the physical Reynolds 
numbers in the systems that one aims to study, but just 
large enough to excite the turbulent small-scale dynamo. 
A strong resolution dependence of turbulent and mag- 
netic properties is thus expected, with higher resolution 
resulting in smaller numerical viscosity and diffusivity, 
higher Reynolds numbers, and thus yielding faster mag- 
netic field amplification. To demonstrate this effect, we 
perform five numerical simulations, in which we resolve 
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the Jeans length, Aj with 8, 16, 32, 64, and 128 cells. 
In our highest resolution simulation with 128 cells per 
Jeans length, we use AMR with 19 levels of refinement 
by a factor of two. This results in an effective resolu- 
tion of 524 288'^ grid cells by the time when we have to 
stop the simulation, because it become computationally 
prohibitive. 

2.3. Analysis in the collapsing frame of reference 

To understand the behavior of the system quantita- 
tively, we need to follow its dynamical contraction in an 
appropriate frame of reference. First, we note that the 
physical time scale becomes progressively shorter during 
the collapse. We therefore define a dimensionless time 
coordinate r (see paper 1), 



1 



tsit) 



dt. 



(2) 



which is normalized in terms of the local free-fall time. 



32G(p(t)) 



1/2 



(3) 



where {p{t)) is the mean density in the central Jeans 
volume, Vj = 47r(Aj/2)^/3. If not otherwise stated, we 
obtain all dynamical quantities of interest within this 
contracting Jeans volume, which is centered on the posi- 
tion of the maximum density. This approach enables us 
to study the turbulence and magnetic field amplification 
in the collapsing frame of reference. However, we also 
analyze radial profiles and magnetic field spectra in the 
fixed frame of reference to study the global evolution of 
the system. 

2.4. Fourier analysis 

To study the scale-dependence of the magnetic field, 
we use Fourier analysis and compute magnetic energy 
spectra, P{B), defined as 



B • B* Ank^ dk , 



(4) 



where B denotes the Fourier transform of the magnetic 
field. The spectral energy density is averaged over spher- 
ical shells in Fourier space, such that the directional in- 
formation of the wavevector is integrated out, and the 
pure scale-dependence of the field can be studied as a 
function of the norm of the wavevector, k = |k|. Note 
that integrating P{B) over all wavenumbers yields the 
total magnetic energy. Thus, P{B) dk is the magnetic 
energy on scales between k and k -\- dk. 

For computing the spectra inside the Jeans volume of 
the core, we extract the AMR data inside a cube with 
length equal to the local Jeans length of the core and in- 
terpolate all relevant quantities for the analysis (density, 
velocity, and magnetic field components) to a uniform 
grid with resolution equivalent to the local Jeans resolu- 
tion. 

A problem in applying Fourier analysis to non-periodic 
datasets is that the spectrum can be distorted by the 
fact that the data is discontinuous at the boundary of 
the extracted cube. Methods to avoid this include ap- 
plying a spherical window function, which smoothly ap- 
proaches zero on all boundaries of the cube, zero-padding 



around the boundaries, or mirroring the dataset on all 
sides. We compared spectra obtained by simply ignor- 
ing the non-periodicity of the data with spectra obtained 
by windowing with a sp herical Hann window, by zero- 
padding (e.g., Brunt et al. 2003), or by mirroring the 
dataset (e.g., |Ossenkopt et ai. 2008a) . We find that the 
spectra are very similar in all cases. They are mostly 
affected at small scales (high-fc) when different meth- 
ods are applied. Since this high-wavenumber range is 
also strongly affected by numerical resolution, we sim- 
ply ignore the non-periodicity of the data and compute 
spectra as if the extracted data were periodic. We also 
tried a wavelet trans f ormation, the A-yariance t echnique 



( [Stutzki et al.||1998[ [Ossenkopf et al.||2008a|b[) to com- 
pute the scale dependence of the magnetic field, and 
found consistent results. However, since the direct spec- 
tral method allows a decomposition of the turbulent ve- 
locity field into solenoidal and compressible modes, and 
thus enables us to study these components separately, we 
prefer the Fourier analysis for the present purpose. 

3. MAGNETIC FIELD AMPLIFICATION BY 
GRAVITY-DRIVEN TURBULENCE 

3.1. Time evolution 

The time evolution of the density, magnetic field 
strength and Mach number of the magnetized, turbu- 
lent collapse of our Bonnor-Ebert sphere is discussed in 
detail in paper I. In particular, we showed that the root- 
mean-squared (rms) magnetic field increases rapidly from 
a small seed field of 10"^ G to the Milli-Gauss level, i.e., 
over six orders of magnitude in our highest resolution 
run with 128 cells per Jeans length. We showed that 
two orders of magnitude of this growth are due to twist- 
ing and folding of magnetic field lines, i.e., due to the 



action of the turbulent sm all-scale dynamo (Branden- 
burg & Subramanian 2005). This is, however, a lower 



limit to the physical amplification by the dynamo, as the 
dynamo growth rate depends strongly on the Rey nolds 
number, and th us on numerical resolution (see, e.g., [Hau- 
gen et al. 2004). If the Jeans length is resolved with 16 
cells or less (which is quite common in typical present- 
day grid simulations involving gravity, and correspond- 
ingly in simulations with th e smoothed particle hy dro- 
dynamics technique, see, e.g., Federrath et al.|2010a ), no 
dynamo amplification is seen at all. Even m our highest 
resolution run, where we resolve the Jeans length with 
128 cells, most of the amplification is due to compression 
of magnetic field lines, i.e., fiux-frcezing, and not due 
to the dynamo. The rather small turbulent Reynolds 
numbers achievable in modern simulations (of the order 
of a few hundred), do not allow for a more efficient dy- 
namo amplification in the simulations. However, since 
the Reynolds numbers in real astrophysical systems are 
typically much higher (e.g., of the order of lO'' in the in- 
terstellar medium) , dynamo amplification will eventually 
dominate over compressional amplification. 

3.2. Radial profiles 

Figure [T[ shows the radial profiles of gas density, rms 
magnetic neld, radial infall velocity, Vr and dispersion, 
<Jv,r, polar and azimuthal velocity dispersions, a^fi and 
<Tv.it> at r = 8, 10, and 12 in the collapse regime for om 
run with 128 cells per Jeans length (note the definition of 
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Fig. 1. — a) Radial density profile of the collapsing, magnetized, 
turbulent Bonnor-Ebert sphere at times r = 8, 10, and 12. b) 
Same as a) but for the rms magnetic field, c) Same as a) but 
for the infall velocity, —v^ and the radial velocity dispersion, (t„ ,.■ 
d) and e) Same as a) but for the polar and azimuthal velocity 
dispersions, u^g and (t„ .0. The vertical lines extending through all 
panels indicate the Jeans radius at t = 8, 10, and 12, respectively. 
See paper I for the corresponding radial profiles of the rms velocity 
and Mach number. 



the diniensionless time coordinate r in eq. [2j) . The den- 
sity and radial velocity profiles follow the t ypical collapse 
profiles of an unstable gas cloud (see, e.g., |Larson||1969[ 
Penston 1969 1. The initial density profile exhibits a short 
power la w with exponent -2.2 for large radii ( |Ebert,1955; 



Bonnor 1956), while the profiles at later times are very 
well fit with a power law, p oc r~^'^. This steepening of 
the profile during the collapse is a result of the slight de- 
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Fig. 2. — Power-law part of the radial profiles of the rms magnetic 
field at T = 12 for Jeans resolutions of 8, 16, 32, 64, and 128 cells. 
Power laws corresponding to pure flux-freezing, Brms c< r"*/"^, and 
a steeper, ex r"^" profile are drawn for comparison. 



viations from an isothermal equation of state (effective 
r = 1.1 instead of 1). The turbulent velocity fluctu- 
ations, measured in terms of the dispersions, av,r, '^v,o, 
and cr„^0 remain subsonic inside the central Jeans volume 
for T < 12, because the sound speed increases during the 
collapse. The rms magnetic field profile exhibits some 
similarities to the density profile, however, the power-law 
exponent in the envelope is different. For the 128 cell run, 
Brias r^^ *^ outside the central Jeans volume, which 
is significantly steeper than what is expected for purely 
compressional amplification, i.e., flux- freezing during the 
collapse (firms oc r^^/^). The steepening is thus due to 
dynamo amplification. However, since the turbulent dy- 
namo is highly resolution-dependent, the radial power- 
law exponent is expected to depend on resolution, which 
we discuss next. 

In Figure [2] we show the power-law part of the radial 
profiles of the rms magnetic field, firms for Jeans resolu- 
tions of 8, 16, 32, 64, and 128 cells at r = 12. Both the 8 
and 16 cell runs are consistent with the radial exponent 
expected for pure flux-freezing, firms r~^^^. Only for 
the runs with 32, 64, and 128 cells, we see a steepening of 
the profiles, which is caused by additional amplification 
of the field through the turbulent dynamo, again showing 
that the Jeans length must be resolved with significantly 
more than 16 cells (e.g., 32 grid cells) for minimum dy- 
namo action to occur. The radial exponent for the 128 
cell runs is close to -2.0, however, the profile is expected 
to steepen further with increasing Jeans resolution. The 
resolution dependence of our results is discussed in detail 
in section [4j where we present evidence for the require- 
ment of a new Jeans resolution criterion in simulations 
of self-gravitating gas. 

Combining the radial dependence of the density, p oc 
and magnetic field, firms oc r"^" in our 128 cell 
run, we find firms oc p^'^^, which is expected to increase 
further at higher Reynolds numbers (higher resolution). 
This is clearly steeper t han the pure flux-freezi ng case. 



fir,: 



oc p 



2/3 



Moreover, 



Schleicher et al. 



'rms OC p'^'-'^ can 



(20091 showed 



ead to a signif- 



that already a scaling of fir 
icant change in the thermodynamics of primordial gas, 
which we show here will have an even stronger effect, 
when turbulent dynamo amplification is taken into ac- 
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count. 

The physical exponent of the radial distribution of the 
magnetic field is set by the physical viscosity and mag- 
netic diffusivity. Since the physical viscosity and diffu- 
sivity can be very small, resulting in Reynolds numbers 
orders of magnitude higher than what we can model in 
a computer simulation, the radial profile of the mag- 
netic field is expected to be significantly steeper than 
Bims oc r^^ ". This means that the potential influence of 
the magnetic field can be significant not only inside the 
Jeans volume, but also on scales larger than the local 
Jeans length, because of the strong amplification of the 
field when the core went through the previous dynamo 
amplification on scales outside the Jeans volume. 

If the Reynolds numbers are sufficiently high inside the 
Jeans core, we would expect that the magnetic field can 
increase to about 10% of equipartition with the turbu- 
lent energy on time scales much shorter than the free-fall 
timcscale, boosting the magnetic field to a dynamically 
significant level even before the core has had time to 
contract much. In the saturated phase, the radial depen- 
dence of the magnetic field is expected to change. The 
saturation behavior, however, cannot be addressed with 
the current simulations, because even at the last available 
time step (after which the simulation becomes computa- 
tionally too expensive to advance any further), at t = 12 
in the 128 cell run, the ratio of magnetic to kinetic en- 
ergy in the core is i^magZ-Ekin ~ 2.2 x 10~^, still far away 
from the expected saturation level (-EmagZ-Ekin ~ 10%). 
Thus, the saturation behavior needs to be investigated 
in a follow-up study with initial field strengths closer to 
equipartition. 

3.3. Turbulence and magnetic field morphology 

To visualize the structure of the density, magnetic field 
and turbulence, we show three-dimensional renderings of 
the density, velocity, magnetic field lines (volume-filling 
and individual ones), the vorticity, |V x v|, and the di- 
vergence of the velocity, V • v, in Figure [Sj The central 
Jeans volume is shown at r = 12, i.e., far into the collapse 
regime, when -through dynamo action and compression- 
the magnetic field has increased by six orders of magni- 
tude to approximately 1 mG. Density fluctuations inside 
the Jeans volume (upper left panel) are rather weak (see 
also the radial profiles in Fig. [T]) and are mostly cor- 
related with the divergence of the velocity field (lower 
right). Since the surroundings of the core are subject to 
contraction, the velocity field is dominated by regions of 
compression. Only some patchy, cloud-like regions show 
expansion. The velocity field (upper right) is turbulent, 
with the outer regions showing some traces of infiow to- 
ward the center of the Jeans volume (see also the ra- 
dial velocity profiles of Fig. [T]). The vorticity contours 
(lower left) are elongated filaments. Some of these fila- 
ments seem to extend through the whole Jeans volume. 
They are folded and twisted several times, similar to the 
magnetic field lines (middle panels). The magnetic field 
is extremely tangled, which is typical of the small-scale 
dynamo. Since the magnetic field is amplified most ef- 
ficiently by winding-up and folding of the field lines on 
scales of the smallest resolvable vortices, the field is most 
strongly twisted on these smallest scales. In the next sec- 
tion, we investigate the scale-dependence of the magnetic 
field quantitatively by means of magnetic and turbulent 



energy spectra. 

3.4. Magnetic field spectra inside the core 

In a proof-of-concept study we suggested in paper I 
that a small seed magnetic field can be amplified sig- 
nificantly by the small-scale dynamo during the collapse 
of a dense gas cloud. Here, we provide a more quantita- 
tive analysis by investigating Fourier spectra to study the 
scale-dependence of the turbulence and magnetic field. 
First, we restrict the analysis to Fourier spectra com- 
puted in the collapsing frame of reference. We therefore 
extract a Cartesian box centered on the maximum den- 
sity in the core with side length equa l to the Jeans length 
at each time frame (see section 2.4 1. 

Figure |4] shows the magnetic energy spectra as a func- 
tion of wavenumber k/kj, i.e., normalized to the local 
Jeans wavenumber, fcj, for our highest-resolution run 
(128 cells per Jeans length). We actually used a three 
times larger extraction volume to include at least some 
information from outside the core, such that all spectra 
were computed on a uniform grid with 384"^ grid cells. 
We did not extract a larger volume, because the data are 
given on an adaptively refined grid, which only contains 
the fiuid variables resolved to the highest level inside the 
Jeans volume. Interpolating the simulation data to sig- 
nificantly larger grids would introduce artifacts. The cho- 
sen volume for extraction onto the highest resolution is 
thus a reasonable compromise to obtain reliable spatial 
information for both the interior of the core and its clos- 
est surroundings. 

An important consistency check for the extraction pro- 
cedure and spectral method is to investigate the initial 
conditions first. Figure |4] includes the magnetic energy 
spectrum at r = 0, which shows that our initial power- 
law dependence, P{B) cx k^^ is reproduced on scales be- 
low the peak. Slight deviations from this power law are 
due to the extraction procedure and Fourier analysis of a 
non-periodic volume (see section 2.4). Strictly speaking, 
the discrete Fourier analysis implies a periodic dataset, 
however, the deviations introduced by simply ignoring 
this are negligible for the present dataset. In addition to 
the power-law scaling, also the peak position must be re- 
produced, which we set to 0.8 pc in the initial conditions 
(see paper I). With the initial Jeans length of 1.15 pc, 
this leads to an expected peak position oik/kj — 1.4, in 
very good agreement with the peak position measured in 
the spectrum at r = 0. 

The peak of the initial magnetic spectrum at about 
k/kj = 1.4 quickly shifts to smaller scales and stays 
roughly constant at fcZ^j ~ 4-6 for t > 4. The initial 
shift of the peak within r < 2 means that the magnetic 
field grows faster on smaller scales, as expected from the 
small-scale dynamo theory. However, the peak does not 
shift further to smaller scales than k/k^ 4-6, which 
corresponds to about 21-32 grid cells, because vortices, 
which drive the dynamo, are under-resolved on scales 



smaller than 30 grid cells (Federrath et al. 2010b). The 
resolution dependence is discussed in detail m section |4j 
On scales smaller than the Jeans scale and larger than 
the peak, the magnetic field spectra are consistent with 
the theoretical prediction of a fc'^/^-Kazantsev spectr um 



( Kazantsev|1968 Brandenburg fc Subramanian|2005 ) in 

the kinematic regime (i.e., the regime in which the mag- 
netic energy is much smaller than the kinetic energy). 
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Fig. 3. — a) Spherical slice of the gas density inside the Jeans volume at t = 12 for our run with 128 cells per Jeans length, b) Velocity 
streamlines on a linear color scale ranging from dark blue (Okms~^) to light gray (5kms~^). c) Magnetic field lines, showing a highly 
tangled and twisted magnetic field structure typical of the small-scale dynamo; yellow: 0.5 mG, red: ImG. d) Four randomly chosen, 
individual field lines. The green one, in particular, is extremely tangled close to the center of the Jeans v olume, e) Co ntours of the vorticity 
modulus, |V X v|, showing elongated, filamentary structures typically seen in subsonic turbulence (e.g., [FVisch|1995| . f) Spherical slice of 
the divergence of the velocity field, V ■ v; white: compression, red: expansion. 
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Fig. 4. — Time evolution of the magnetic field spectra for the 
run with 128 cells per Jeans length as a function of wavenumber, 
normalized to the local Jeans wavenumber, i.e., the magnetic field 
inside the dense core is shown for each timestep, in the collapsing 
frame of reference. The initial magnetic field spectrum (r = 0) 
follows the initial power law, P{B) oc on scales smaller than 
the peak scale at fc/fcj 1.4. Within the regime of turbulent decay 
(0 < r < 4), the peak shifts to k/kj ^ 5-6 due to dynamo ampli- 
fication, which is most efficient on small scales (discussed further 
in section [4. 1[ |. The magnetic field spectra are consistent with the 

faerie k^^'' Kazantsev spectrum I Kazantsev|[l968, BrandenburgI 
Subramanian 2005 ) on scales smaller tJian trie Jeans length and 
ger than the peak. Within the collapse regime, the magnetic 
energy grows exponentially on all scales within the core. 



Indications of the Kazantsev spectrum are also foun d in 
simulations of the intra-galaxy cluster medium by |Xu| 
et al. (2009 2010). However, in both their simulations 
and in ours, the scaling range for measuring the slope 
is too narrow to draw final conclusions. It is not clear 
whether the A;^/^-power law would persist at higher reso- 
lution. In turbulence-in-a- box simulations with external 

(|2004|) with up to 512^ grid 



forcing by Haugen et al 

cells, there is strong indication that the magnetic en 
ergy spectra indeed converge to the Kazantsev slope with 
increasing resolution. As in previous turbulence-in-bo x 
calculations with external forcing (e.g., Cho et ar]|2009 ), 
we find that the magnetic energy grows exponentially on 
all scales, which is a typical feature of the small-scale 
dynamo (Brandenburg & Subramanian 2005), but has 
(to the best of our knowledge) not been shown before 
in a gravity-driven turbulent gas core. We also measure 
the exponential growth rates of the magnetic field as a 
function of Jeans resolution in section |4] below. 



3.5. Magnetic field spectra in the fixed frame of 
reference 

The magnetic energy spectra in Figure [4] suggest that 
the spectra fall off more steeply than the k^^^ Kazant- 
sev law toward lar ge scales, outsid e the Jeans volume (see 



also the spectra in Xu et al |2010 ). In order to clarify this 
and to investigate the field structure outside the Jeans 
volume, we introduce and apply a new method to infer 
the spectrum over more than four orders of magnitude 
in length scales. Since we aim to investigate the field 
structure in the fixed frame of reference at the largest 
available scales in the simulation, but at the same time 
would like to include spectral information on the very 
smallest scales, we apply a two-step approach. First, we 
re-normalize the high-resolution spectra obtained inside 
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Fig. 5. — Magnetic energy spectra in the fixed frame of reference, 
i.e., the spectra are normalized to the initial Jeans wavenumber, 
fcjo . The time evolution in the collapse regime shows that the mag- 
netic field is basically frozen-in on large scales, while the turbulent 
dynamo keeps amplifying the field on smaller and smaller scales as 
the collapse proceeds. 

the Jeans volume (Fig. [4]) , by shifting them to the cor- 
rect position with respect to the fixed frame of reference. 
However, the spectra at late times during the collapse, 
obtained with this method, do not contain any informa- 
tion on scales far outside the Jeans volume. Thus, in 
the second step, we add spectral information on large 
scales, by gradually extracting larger boxes, centered on 
the core at a fixed grid resolution. Using this method, 
we can test the spectral energy scale-by-scale. We call 
this method 'scale-by-scale extraction'. 

The result of this two-step approach is shown in Fig- 
ure [5] where we plot the spectra of magnetic energy in 
the fixed frame of reference, i.e., as a function of the 
initial (at r = 0) Jeans wavenumber, fcjo, of the col- 
lapsing system. The first thing to note is that the spec- 
tra on large scales, obtained with the scale-by-scale ex- 
traction method connect reasonably well with the spec- 
tra computed inside the Jeans volume and re-normalized 
to match the fixed frame of reference. Some devia- 
tion is seen only at k/kjQ sa 300 for the spectrum at 
r — 12, which can be taken as a measure of the un- 
certainty in the spectra obtained with our scale-by-scale 
extraction method. Given the total range of scales that 
we aim to probe here, the difference in the spectra ob- 
tained with our two-step approach is acceptable. We 
also tested whether extending the scale-by-scale extrac- 
tion to smaller scales (inside the Jeans volume) matches 
the high-resolution spectra of step one. We found that 
they do within an uncertainty of about 25%, i.e., the 
slopes and peak positions are reproduced reasonably well 
with the scale-by-scale extraction method. However, the 
high-resolution spectra inside the Jeans volume are more 
accurate, and we thus prefer the two-step approach de- 
scribed above. 

Three main results can be extracted from Figure [51 
First, the magnetic field always grows fastest on the 
smallest available scales in the simulation, with the peak 
located at around 30 grid cells. Second, as the collapse 
proceeds, the field is amplified on small scales inside the 
core, but is essentially frozen-in on large scales. Third, 
the large-scale spectra follow a power law close to fc'^ ". 
This particular exponent for the power law is consistent 
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with the radial dependence of the field (see Figs. [T] and |2| 
and depends strongly on resolution. The magnetic en- 
ergy on scales outside the core is 



cx 



J P{B) dk cx /c^ " . 



(5) 



From this it follows that i?rms « ° oc r~^ '^ outside the 
Jeans volume, which is consistent with the power-law be- 
havior of the radial profile of i?rms in Figure [l] It should 
be re-emphasized, however, that the exact exponent of 
this power law is essentially meaningless, as it depends 
on the numerical Jeans resolution. Higher resolution will 
lead to a steeper increase of P{B) toward small scales 
(see Fig. [2]). In reality, the magnetic field will grow much 
more strongly due to dynamo action than what we can 
resolve in the present calculation with 128 cells per Jeans 
length, and thus, our amplification rate is a lower limit 
(discussed further in section |4|. 

3.6. Probability distribution Junction of the gas density 

The probability distribution function (PDF) of the gas 
density is a useful measure of the turbulence in any tur- 
bulent system exhibiting density fluctuations. Moreover, 
the PDF is an essential ingredient for models of star for- 



mation (e.g.,'Padoan & Nordlund'2002','Krumholz & Mc- 



Kee 2005 



Henncbellc & Chabricr 2008, 2009) and the gas 
distribution in galaxies (e.g., |Tassis|20 07 Krumholz et al. 
|2009J . The density PDF has been studied in some detail 
m non-self-grav i tating, turbulent systems (e.g., Vazquez- 
Semadeni||1994[ [Padoan et al.||1997[ [Passot &: Vazquez- 



Semadenih998 



1 19941 [F 
19981 |L 



emaster &: Stone|2008 Federrath et a. 



2008E' [Price fc Federrath'2010'), and in turbul ent systerns 
includ ing self-gravity (e.g., Klessen 2000; Federrath et al. 



2p08a; Kainulaincn et al.|20D9|iChcrfc Kim|20lH|K.ritsug 

However, it has not been analyzed yet in a 



\et al.j2011 



collapsing system, in which turbulence is replenished by 
the gravitational collapse of a dense core. 

Figure [6] shows the time evolution of the PDFs of the 
logarithmic density 



In 



{P) 



(6) 



where {p) denotes the mean density in the core, i.e., in- 
side the Jeans volume. The PDF at t = is purely 
a result of the initial, radial density distribution, fol- 
lowing a Bonnor-Ebert profile. This profile exhibits a 
flat inner core for radii smaller than the Jeans radius, 
and can be approximated with a power l aw of the form 
p (x r ~° with a ~ 2.2 at large radii ( Ebcrt|l955 |Bon nor| 
[1956 ). Using the derived relation between a power-law 
radial distribution and the corresponding density PDF 
(see Appendix |Aj) , we can estimate the power-law expo- 
nent of the density PDF from the power-law exponent 
of the radial distribution, a. Thus, the PDF of the log- 
arithmic density, s, follows a power law for small loga- 
rithmic densities with exponent — 3/a, and falls off more 
steeply toward higher densities, due to the flattening of 
the Bonnor-Ebert profile in the center of the core (see 

Fig.m). 

Both in the regime of turbulent decay (t < 4), and in 
the collapse regime (r > 4), the volume- weighted PDF 
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Fig. 6. — Time evolution of the probability distribution functions 
(PDFs) of the logarithmic gas density, s = ln(p/ (p)), inside the 
collapsing core. The initial density PDF is characterized by the 
asymptotic radial profile of a Bonnor-Ebert sphere at low densities, 
p (X r~^'^, while at higher densities, it turns into a steeper PDF due 
to the flat inner core. The PDFs in the turbulent decay regime (r < 
4) are convolutions of the PDF of the initial density distribution 
and log-normal distributions typical of turbulence. In the collapse 
regime (t > 4), the PDFs are close to log-normal distributions, 
with some modification due to the radial density distribution (see 
Fig.n). The dotted curve shows a log-normal fit using equation ([tJ 
with rhe imposed condition, eq. Jsl at r = 10. 



develops a log-normal form, 
1 



Ps{s)ds = 



exp 



.turb 



2a2 



;,turb 



ds , 



(7) 



where (s) and Us.tmh denote the mean logarithmic den- 
sity and standard deviation, respectively. This PDF 
is typical of compressible, nearly isothermal, turbulent 
fiows, wh ich has been motiv ated with the central limit 
theorem ( |Vazquez-Semadeni||1994) . Note that the mean 
is related to the standard deviation by 



1 



s.turb ' 



(8) 



because the mean density inside the Jeans volume must 
be recovered by integrat ion of the PDF, weighted bv gas 



density, (p) = f p Psds ( Vazquez-Semadeni 



1994 



Feder- 



rat h et al.||2008b[ |2010b[ ). I'he measured PDI-'s m I'lg- 
ure [Sj however, exhibit significant deviation s from the 



log-normal function due to i ntermittency (.Kritsuk et al. 
2007| [Schmidt et al. 2009 Fe derrath et al.[2010b[ ) and the' 
convolution with the radial density distribution inside 
the Jeans radius of the collapsing Bonnor-Ebert sphere. 
Note that we do not see the typical power-law PDFs 



observed in s elf-grayitatin g systcnis (Klessen 2000j Fed 
errath et a l. 2008a; 'Kai nulainen et al . 2009; C ho &: Kim 
2011| I^riTsuk et al.,,2011 ) during the collapse, because 



here we always compute PDFs in the collapsing frame of 
reference, thus ignoring the gas outside the Jeans volume, 
which has a clear power-law PDF. 

In order to compare the PDFs obtained in our collaps- 
ing system with results on the PDF reported in the liter- 
ature, we compute the standard deviations of the PDFs 
below. The standard deviations of the density PDF are 
as = 0.225, 0.429, 0.353, and 0.290 for r = o', 1, 2, and 
3 in the turbulent decay regime. It is important to note 
that the initial dispersion at r = has nothing to do with 
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turbulence. It is purely due to the initial, radial density 
distribution of the Bonnor-Ebert sphere. In order to sep- 
arate the standard deviation caused by turbulent com- 
pression from the standard deviation due to the radial 
density distribution inside the Jeans volume, we subtract 
the initial density dispersion caused by the Bonnor-Ebert 
profile: <J1^^^^y^ — ctUt) — crs(O). The standard devia- 
tions of the turbulence, Os turb, can then be compared to 



jPadoan et al. 


1997; Passot & Vazqucz-Scmadeni 


1998 


|Ffederrath et al 


. 2008b, 2010b, Price ct al. 2011 

11 ! , i 1 11 


): 



s,turb 



ln(l -h6^A^^ 



(9) 



The parameter h in this equation is of order unit y and de- 
pends on t he way of turbulent energy injection (Schmidt 
et al.|2009|). It can vary by a factor of three, froni b « 1/3 



tor turbulence excited by solenoidal (rotational) modes 
to & « 1 for turbulence^^ excited by purely compressive 



modes (Federrath et al 
smoothly from 1/3 to 



2008b I. Moreover, h mcreases 
, when the mixture of modes 



for the turbulence driving is s moothly varied from fully 
solenoidal to fully compressive (Federrath et al. 2010b). 
Given the rms Mach numbers ot order unity (see paper 
I), we obtain h w 0.42, 0.40, and 0.33 for r = 1, 2, and 
3, respectively, with an uncertainty of ±0.05 each. The 
values at early times (r = 1, and 2) fit the expected value 
(6 K, 0.4) for a natur al mixture -which were our initial 
conditions- very well ( [Federrath et al.||2010b| ). The fur- 
ther evolution indicates a transition to a more solenoidal 
behavior of the turbulence. This is reasonable, as 
decreases from unity to about 0.5-0.6 at t = 3, which 
means that shocks are basically absent at this stage. 
Some caution should be exercised in deriving h for the 
subsonic regime, although equation ([9]) seems to hold 



in both the sub sonic and transsonic regimes (Passot & 
Vazquez-Semadcni„1998| , a s well as in the hi ghly super- 



sonic regime ot turtulence ( Price et"ar]|2011 ). 

In the collapse regime (r ^4), it becomes increasingly 
difficult to separate the density dispersion caused by the 
infall profile and caused by the turbulence that is gen- 
erated during the collapse. The density PDFs are close 
to log-normal distributions in the collapse regime, with 
some modifications due to the infall profile. As an ex- 
ample, we show a log-normal fit to the data at r = 10 
in Figure Ig] as the dotted line labeled p(s). Using the 
total standard deviations, CTs, and the Mach numbers in 
the collapse regime, we derive h « 0.65, 0.48, and 0.37 for 
r = 8, 10, and 12, respectively, again with an uncertainty 
of about ±0.05. These values are upper limits, because 
of the additional dispersion from the radial density dis- 
tribution, which is, however, quite small inside the Jeans 
volume (see Fig. [I]) . These results for the collapse regime 
are indicative of a time evolution of h from a state, when 
compressive modes are more important (at r sa 8) to 
values close to equipartition at later times, which is con- 
firmed with a spectral decomposition into solenoidal and 
compressible modes of the velocity field in section |4.2| 
below. 

3.7. The driving scale and spectrum of gravity-driven 
turbulence 

As discussed in the introduction, the driving of tur- 
bulence by extracting potential energy from the gravita- 



tional field o f a gas cloud has been suggested in several 
studies since Hoyle (1953). But what is the characteris- 
tic driving scale ot gravity-driven turbulence? Answering 
this question is important, because the injection scale of 
turbulence may determine the mode of star formation 
in molecul ar clouds, isolated versus clustered (Klessen] 
et aH2000[[Heitsch et al.|2001[ [Klessen 200l|. Moreover, 



kinetic energy injected on large scales can produce large- 
scale coherent gas compressions, promoting the forma- 
tion of condensations, while turbulence driven on small 
scales acts to inhibit star formation. Over the last 30 
years it has become clear that the turbulence in present- 
day molecular clouds contains mo st of the energy on large 
scales (Larson 1981; Solomon et al. 1987 Ossenkopf &:[ 
Mac Low|2002i iHeyer fc Brunt I2DD11 ), i.e., it is driven on 



large scales (Brunt et ai. 2U09[ an? references therein) 



Most important drivers oi galactic-scale turbulence are 
the expansion waves produced by supernova explosions 
and expanding HII regions around massive star cl usters 
(jMac Low & Klessen 120041 [Schneider et al.||2010|), and 



the galac tic spiral shock^ producing large-scale converg- 



ing flows dVisim iac 1994'; Wa lder fc Fo lini 1996; Vazquez- 



Semadeni et al. 2006; Hcitsch et al. 2006 ; Hcnnebclle 
et al 112008; Ba nerjeeet al. 2009 ; Folini'eral. ||20lU| . The 



power of jets and outflows from low-mass stars is also 
enormous, and may thus contribute to the tur bulence 
driving on intermediate and small scales 



Tan & 



iBlackman 2004; Nakam ura fc L i 2008; Li et al.||i^OTnf 
Carroll ct al. 2010). 

'lb investigate the driving scale and scale dependence 
of our gravity-driven turbulent core, we compute veloc- 
ity spectra inside the Jeans volume of our 128 cell run, 
defin ed in analog y to the magnetic field spectra above 
(e.g., |Frisch||1995| as 



P{v) dk ■ 



V • V* 47rfc^ dk , 



(10) 



where v denotes the Fourier transform of the velocity 
field. The total Fourier spectrum can be separated into 
transverse (k ± v) and longitudinal (k || v) parts by ap- 
plying a Helmholtz decomposition. Note that integrating 
the transverse velocity spectrum yields the specific ki- 
netic energy in transverse (rotational) modes, while inte- 
grating the longitudinal velocity spectrum yields the spe- 
cific kinetic energy in longitudinal (compres sible ) modes. 
We will use this decomposition in section |4.2[ to inves- 
tigate the relative mixture of modes in the turbulence 
spectrum. First however, we concentrate on the total 
(transverse + longitudinal) velocity spectrum. 

Figure [Tl sh ows the velocity Fourier spectra, defined in 
equation p_0| ). In the upper panel we show the spectra 
for the fulTvelocity field, while the bottom panel shows 
the spectra, where we subtracted the radial infall profile 
(see Fig. fTl) from the velocity field before computing the 
spectra, this way, we obtain the pure spectrum of 
turbulence without direct contributions from infall. As 
for the magnetic field spectra, the initial velocity power 
spectrum is reproduced with our Fourier analysis. While 
the turbulence is decaying (t < 4), the initial spectrum 
following f lattens slightly and turns into the Kol- 

consistent with 



mogorov 



( 1941 ) spectrum, P{k) cx k" 
the expectation for subsonic, hydrodynamic turbulence 
(e.g., )Frisch|,1995| ). For MHD turbulence, a similar seal- 
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Fig. 7. — Top: Velocity power spectra normalized to the local 
Jeans wavenumber, k/kj, inside the core for different times until 
T = 12. Bottom: same as top, but for the infall-subtracted veloc- 
ity spectra. The flattening of the infall-subtracted spectra around 
fc/fcj ^ 1-2 indicates that gravity-driven turbulence exhibits an 
effective driving scale close to the Jeans scale. 



ing exponen t is found (e.g., Cho & Vishniac 2000 Cho 



et al.||2009 [). Since the turbulence inside tfie Jeans vol 
ume remains subsonic for r < 12 (see paper I), the Kol- 
mogorov scaling persists until the end of the simulation. 
The scaling range, however, only extends over the in- 
terval 2 < fc/fcj < 5, because of the rather low resolu- 
tion inside the Jeans volume (128 cells). Measuring the 
power-law exponent of the turbulence accurately would 
require at least 512 cells inside the Jeans volume as ex- 
pect ed^ran_2ure_turbu^^ 



2010b 



Price fc Federrath||2010 |. 



box (Federrath et al 

In the collapse regime, the spectra retain their over 
all shape, but shift upward with increasing time, which 
means that the specific turbulent energy, and thus the 
velocity dispersion increase in the collapse regime (see 
Fig. 1). The infall-subtracted spectra (bottom panel of 
Fig. 7 1 become quite flat around the Jeans scale, which 
indicates that the turbulence inside the core is driven 
from the outside with an effective driving scale, ^int, 
approximately on the Jeans scale, Aj. Since the Jeans 
scale continuously decreases during the collapse, gravity- 
driven turbulence does not have a fixed driving sc ale 



This confirms the assumption in Schleicher et al. ( 2010 1 
that the integral scale of the turbulence is a traction of 
the Jeans scale in their model of dynamo-driven magnetic 
field amplification during the formation of the first stars 
and galaxies. They assumed £int = O.lAj, while we find 



05 

a:" 
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Fig. 8. — Shows the resolution dependence of the magnetic field 
spectra at r = 12. The peak of the magnetic energy shifts to 
smaller scales for increasing Jeans resolution (32, 64, and 128 cells), 
but is completely absent for a Jeans resolution of 16 cells and below. 
We conclude that turbulent dynamo amplification of the magnetic 
field is only obtained, if the Jeans length is resolved with at least 
30 grid cells. 

in the present numerical experiment that the effective 
turbulence driving scale is somewhat closer to the Jeans 
scale, £int < Aj. 

4. A NEW JEANS RESOLUTION CRITERION FOR 
(M)HD SIMULATIONS OF SELF-GRAVITATING 
GAS 

In this section, we analyze the resolution dependence 
of the turbulence and magnetic field growth inside the 
Jeans volume of our collapsing, magnetized core in more 
detail. We show that a Jeans resolution of 16 cells is 
clearly insufficient to obtain dynamo amplification of the 
magnetic field. Only with a Jeans resolution of 32 cells, 
we see dynamo amplification, which we explain with the 
fact that rotational motions of the turbulence (e.g., vor- 
tices) are not well resolved, if the Jeans length is re- 
solved with 16 cells or less. In contrast, if the Jeans 
length is resolved with 32 cells and more, the turbulent 
energy in the core is converged with resolution. We thus 
conclude that a minimum Jeans resolution of about 30 
grid cells (depending on th e particular numerical scheme, 
see, |Kitsionas et al.||2009|) is required to sufficiently re- 



solve the turbulent energy, and thus the turbulent pres- 
sure, and to obtain minimum turbulent dynamo action 
close to the resolution limit. We emphasize that -apart 
from a few exceptions- almost all present-day simula- 
tions of self-gravitating magnetized and non-magnetized 
gas use a Jeans resolution criteri on significantly below 
our threshold value of 30 cells (e.g., Truelove et al.|1997 ). 



4.1. Resolution dependence of the turbulent dynamo 

The time evolution of the rms magnetic field in pa- 
per I and the radial magnetic field profiles of Figure |2] 
demonstrate a strong resolution dependence of the tur- 
bulent dynamo. In this section, we explain the resolution 
threshold for dynamo action, which we find is about 30 
grid cells per Jeans length (see also paper I). Figure Js] 
shows the resolution dependence of the magnetic fiela 
spectra at t = 12. Our five runs, in which we resolved 
the Jeans length with 8, 16, 32, 64, and 128 grid cells are 
compared. The runs with 32, 64, and 128 cells show a 
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Fig. 9. — Shows the time evolution of the rms magnetic field 
strength as a function of Jeans resolution, where we divided out 
the effect of pure flux freezing, which can provide a maximum pos- 
sible field amplification by p^^^ inside the core. The remaining 
amplification is thus due to the small-scale dynamo. The num- 
ber on each curve indicates the growth rate, Q, in the expression 
Brms/ p^^^ oc exp(nT), measured in the interval t = [8, rend]- 



clear maximum of the magnetic energy density on scales 
below the Jeans scale, i.e., inside the core region, while 
the 8 and 16 cell runs do not exhibit such a peak. The 
peak shifts to smaller scales with increasing Jeans resolu- 
tion: it is located at k/kj « 1-2 for 32 cells, fc/fcj « 2-4 
for 64 cells, and fc/fcj « 4-6 for 128 cells. The peak in 
P{B) thus always appears close to the resolution limit, 
on scales of about 20-30 grid cells, which means that the 
dynamo amplification of the magnetic field is strongest 
on these scales. This result, taken together with the ab- 
sence of a peak for the run with Jeans resolution of 8 
and 16 cells suggests that a Jeans resolution of at least 
30 grid cells is required for the dynamo to work in a grid 
simulation. This is significantly more than the resolution 
criterion of 4 grid cells per Jeans length to a void a rtifi- 
cial fragmentation found by Truelove et al. (1997) and 
frequently applied in simulations involving selt-gravity 



with grid codes (e.g., Krumholz et al.||2004|) and corre- 



spondingly with pa r ticle codes (e.g. , Bate et al. 1995; 
[Price fc Batel[20U7{ [Federrath 2010a, and refer- 

ences therein). The explanation for this is that magnetic 
field amplification by the turbulent dynamo is most effi- 
cient on the smallest scales due to the fast eddy turnover 



times on t he smallest scales (e.g., Brandenburg & Sub- 
[ramanianl 2005 , for a review on turbulent dynamo am- 
plification)^ However, the dynamo feeds from the trans- 
verse, solenoidal turbulent motions, which -due to the 
discretization of the fluid variables- are only resolved, if 
they have at least 30 grid cells across. This is consis- 
tent with the spectral analysis of high-resolution (1024'^) 
fi xed-grid simulations of driven, supersonic turbulence 
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(2010b), where a clear deficit of rota- 
tional turbulent eiiergy was detected below 30 grid cells. 

Figure [9] shows the time evolution of the magnetic 
field for oifferent Jeans resolutions. The plot shows 
Brms/ P^^^ , i-e., we divided out the maximum possible 
amplification by pure fiux-freezing in spherical symme- 
try. In this representation, we can measure the field 
growth caused by the turbulent dynamo (see paper I). 
The growth rate, f), in the expression B^^g/p' 



o 




1000 



Reynolds number 



Fig. 10. — Dynamo growth rate, f2, measured in Fig.|9]as a func- 
tion of Reynolds number. Note the strong increase in the growth 
rate for Rm < 40, which corresponds to a Jeans resolution of about 
30 grid cells. 



exp (rir), is indicated as a label on each curve. For 8 and 
16 cells, no clear dynamo amplification occurs, while the 
growth rate increases discontinuously to = 0.29 for a 
Jeans resolution of 32 cells, marking the onset of dynamo 
action. 

We note that the initial decrease of i?ims/p^^'^ for the 
runs with 8 and 16 cells for r ;< 6 is mostly due to the 
turbulent decay and not strongly affected by numerical 
diffusion of the magnetic field. We estimate the effects 
of numerical diffusion in Appendix [C] which may ac- 
count for a deviation from ideal MHD by e < 0.005 in 
B cx p'^/^~'^ for any of the runs. The 8 and 16 cell runs 
show a significantly lower level of turbulence than all the 
other runs (see paper I, Fig. 3e), because the turbulent 
energy is not sufficiently resolved in those runs, which is 
discussed in detail in section [42] 

To see the discontinuous increase of the growth rates 
with threshold resolution of about 30 cells more clearly, 
we plot the growth rates as a function of Reynolds num- 
ber in Figure 10 For computing the Reynolds number, 
we assumed that the most dissipative wavenumber in the 
simulation, kr, = N/2, corresp onding to two grid cells 



(similar to Haugen et al.[[2004 ): 



4/3 



(11) 



with the number of grid cells per Jeans length N. Note 
the str ong increase in the growth rate for Rm < 40 in Fig- 
ure 10 For Rm > 40, the increase in the growth rate can 
be approximated with a power law, c>c Rm"'^, slightly 
shallower than the theoretical dependence on Reynolds 
number (fi cx Rm°'^, see Appendix|B]). In turbulence-in- 
a-box studies without gravity, a critical Reynolds number 
of about 35 was foun d for the onset of dynamo action 
in excellent agreement with the 



( Haugen et al 
present study o: 



2004[ ) 

a gravity-driven turbulent core. 



oc 



4.2. The mixture and resolution dependence of 
turbulent compressible and solenoidal modes 

In Figures [2] [8] and [9] we showed that the small-scale 
turbulent dynamo only starts to operate, if the Jeans 
length is resolved with about 30 grid cells. Why is 
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Fig. 11. — Shows the resolution dependence of the velocity spec- 
tra at r = 12. The spectra do not show convergence, even with a 
Jeans resolution of 128 cells. We expect that a resolution of 512'^ 
grid cells inside the Jeans volume is req uired to obtain an inertial 
range. The expected |Kolmogorov| |jl941 ) scaling of the turbulence 
in the inertial range is shown only tor comparison. The total tur- 
bulent energy, i.e., the integral over the spectra, however, indi cate s 
convergence for a Jeans resolution of 32 cells and more (see Fig.|12|. 



there such a threshold resolution? In this section we 
show that the solenoidal (rotational) turbulent motions, 
which drive the dynamo are severely under-resolved, if 
the Jeans length during the collapse is resolved with 16 
cells or less. On the other hand, the solenoidal turbulent 
energy converges for 32 cells per Jeans length and higher 
resolution, indicating a threshold resolution between 16 
and 32 cells per Jeans length for the dynamo to operate. 

First, we plot the resolution dependence of the infall- 
subtracted velocity spectra (equation [lO| in Figure 11 
These spectra are not converged. Although the 128 



cell ru n seems consistent with the expected |Kolmogorov| 
( 1941 1 scaling for subsonic turbulence on some length 
scales, there is no true inertial scaling range. This is not 
surprising, given that previous turbulence-in-a-box sim- 
ulations of driven turbulence showed that inertial range 



scaling requires at least 512 cells (e.g., Haugen et al.|2004 
Kritsuk et a"L][2007l [Schmidt et aL|[2009; Fcderrath et al'. 
2010bD. Although the spectral form is clearly not con- 



verged at a Jeans resolution of 128 cells, the integral 
over the velocity spectrum, i.e., the turbulent energy or 
turbulent pressure shows indeed convergence for a Jeans 
resolution of 32 cells and higher. 

Figure [12] (top panel) shows the solenoidal ratio in- 
side the gravity-driven turbulent core. We define the 
solenoidal ratio as the specific kinetic energy in rotational 
turbulent motions divided by the total specific kinetic en- 
ergy on scales smaller than the Jeans scale, 



X 



Xt7 Psoi{v)dk 
E^t ~ P{v)dk 



(12) 



The solenoidal ratio at r = is y = 2/3, which is 
the na tural mixture ( Elmegreen fc Scalo|2004{ |Federrath| 
et al. |2008b) , imposed as an initial condition! About 
2U7o 01 the initial compressible modes are dissipated very 
quickly, because the turbulence is decaying, and thus x 
increases from 2/3 to about 0.8 in the turbulent decay 
regime, r < 4. In the first stages of the collapse regime 
(4 T 7), X drops to about 40-50%. Within the 




til 



3 
2.5 

2 
1.5 

1 

0.5 




128 cells 

64 cells - » - 

32 cells - — • 

16 cells ---B--- 

8 cells -■- 

radial infall profile subtracted 




2 4 



10 12 14 16 18 



Fig. 12. — T op: time evolution of the solenoidal ratio, x = 
Eaoll Etot, eq. \12\ , for Jeans resolutions of 8, 16, 32, 64, and 
128 grid cells. Bottom: Turbulent energy in solenoidal modes for 
the infall-subtracted motions inside the core. The runs with 32, 
64, and 128 cells indicate convergence of x (top panel) and in the 
absolute rotational energy, E^^i (bottom panel), while the 8 and 
16 cell runs clearly underestimate the amount of rotational energy. 



collapse phase, the kinetic energy provided by the grav- 
itational collapse is gradually converted into rotational 
motions, which drive dynamo amplification of the mag- 
netic field. For r > 12, the solenoidal ratio seems to 
approach 2/3, which correspon ds to the natural ratio in 



three-dimensional turbulence (Federrath et al. 2008b) 



This natural ratio is most easily pictured by recalling the 
number of spatial directions in which a longitudinal wave 
can compress the gas in a three-dimensional system. It 
is 1 out of 3, which leaves 2/3 for rotational motions. 

Gravity-driven turbulence inside galactic disks i n two- 
dimensional (2D), high-resolution simulations by |Wada 
et al. (2002) show comparable values of x ~ 1/2 to what 
we find here. In 2D, x is indeed expected to be 1/2 from 
the dimensional arguments above. However, 2D turbu- 
lence is different from three-dimensional (3D) turbulence, 
as there is an inverse cascade in 2D, i.e., vortices merge 
to form larger vortices instead of breaking up into smaller 
vortices as in the 3D case. Thus, the turbulent dynamo 
is expected to behave (if at all present) very differently 
from the 3D case studied here. 

The evolution of the solenoidal ratio toward more 
solenoidal turbulence for r > 7 explains that the growth 
rates in Figure ]9] increase after t — 7, showing a steepen- 
ing of i?ims/ /O^'^with time for r 7-11, which correlates 
with the increase in the solenoidal ratio. This is because 
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the dynamo feeds from the solenoidal modes of the tur- 
bulence, which increase between r = 7 and 11. Since 
the solenoidal ratio varies strongly with time in the early 
phases of the collapse, we had to choose a rather late time 
within the collapse regime to measure the growth rates 
of the magnetic field, i.e., we estimated the growth rates 
in the interval t = [8, Tend] in Figure [g] Ideally, we would 
like to measure the growth rates after the solenoidal ratio 
has converged to x = 2/3, but due to the computational 
expense, we could not follow the high resolution simula- 
tions (64 and 128 cells) far into this regime. The 32 cell 
run, however, does show a convergence of x toward 2/3. 

Both the solenoidal ratio and the absolute specific ki- 
neti c en ergy in rotational motions, Egoi, shown in Fig- 
ure [l^l indicate convergence for a Jeans resolution of 
32 cells and higher. Using 8 and 16 cells to resolve the 
Jeans length is clearly insufficient to resolve rotational 
motions. We conclude that at least 30 grid cells per 
Jeans length must be used in (M)HD simulations of self- 
gravitating gas to resolve the turbulent energy, and thus 
the turbulent pressure, and to obtain minimal dynamo 
amplification of the magnetic field on the Jeans scale. 

5. CONCLUSIONS AND IMPLICATIONS 

We presented high-resolution magnetohydrodynamical 
simulations of the collapse of a dense, magnetized gas 
cloud (see also, Sur et al.|2010 paper I) . During the col- 
lapse of the cloud, gravitational, potential energy is con- 
verted into turbulent motions, which in turn amplify the 
magnetic field exponentially fast by the turbulent dy- 
namo process. The exponential amplification is driven 
by the stretching, twisting, and folding of the small- 
scale magnetic field lines (Fig. Isl). At sufficiently high 
Reynolds numbers, even extremely small initial seeds of 
the magnetic field are expected to be amplified to dynam- 
ically important magnetic field strengths on timescales 
much shorter than the collapse timescales. We conclude 
that magnetic fields should not be neglected in both pri- 
mordial and contemporary studies of star formation. 

We investigated the scale-dependence of the turbu- 
lence and magnetic field by means of Fourier analysis in 
the collapsing frame of reference (Fig. |4]) , showing some 
indication of the Kazantsev spectrum of turbulent dy- 
namo amplification, and in the fixed frame of reference 
over more than four orders of magnitude in spatial scale 
(Fig. [5|. Wc find that the effective kinetic energy in- 
jection scale of gravity-driven turbulence is close to the 
Jeans length (Fig. ItI. 

Our Fourier analysis of the magnetic field shows that 
the dynamo is only excited, if the Jeans length is suffi- 
ciently resolved (Fig. [sl . The radial dependence of the 
magnetic field is signmcantly steeper than what is ex- 
pected from purely dragged- in magnetic field lines, i.e., 
flux-freezing (i?rms r^^^'^). For a Jeans resolution of 
128 cells, we obtained Brms c>c r"^'*', which is expected to 
steepen further with increasing Jeans resolution (Fig. [2]) . 
We find that at least 30 grid cells per Jeans length are 
required for minimum dynamo action to occur in col- 
lapse simulations, a resolution requirement, which is not 
achieved in most current simulations. Here, we studied 
dynamo amplification with a resolution up to 128 cells 
per Jeans length, while usually less than 16 cells are used. 
The amplification rate increases with resolution (Fig. [o]) , 



which renders any existing simulation of dynamo ampli- 
fication of a highly turbulent medium insufficiently re- 
solved to determine the physical growth rate of the mag- 
netic field, and can at best provide lower limits on the 
physical growth rates. 

We find that the probability distribution function of 
the gas density inside the gravity-driven, turbulent Jeans 
volume follows a log- normal distribution (Fig. l6]) . The 
standard deviations and Mach numbers inside the core 
indicate that gravity-driven turbulence approaches a 
natural mixture of solenoidal and compressible modes, 
Eso\/Etot ~ 2/3, after a phase of more compressively 
driven turbulence, caused by the global collapse of the 
system (Fig. [T2| . The turbulent energy (or turbulent 
pressure) converges only for a Jeans resolution exceeding 
30 grid cells. In contrast, the solenoidal component of the 
turbulent energy is severely under-estimated, if the Jeans 
length is resolved with 16 cells or less. Thus, we suggest 
a new Jeans resolution criterion of 30 grid cells per Jeans 
length to obtain converged results of the turbulent energy 
on the Jeans scale, and to capture minimum magnetic 
field amplification by the turbulent dynamo process. 

The importance of magnetic fields in present-day ac- 
cretion disks is generally accepted. However, since even 
small initial seeds of the magnetic field are efficiently am- 
plified by turbulent dynamo action, it cannot be excluded 
that magnetic fields also play an important role in pri- 
mordial accretion disks. Simulations of primordia l gas 
show that it is highly turbulent (e.g. , Abel et al. 20021 
O'Shea fc Normaii|[2007l [Wise fc Abel, 2007. ; Clark ct at 



2om 



cient 



et al. pOOSf TwhicF suggests that "tkere is suffi- 



tmetic energy m rotationa l modes of the turbulence 
(e.g., vortic ity; see in particular, [Wise fc Abeil2007} [Crei^ 
et al. 2008) to drive the small-scale dynamo als o In pri 



mordlal gas clouds and accretion disks (see also, Pudritz 
fc Silk||T989| ) 



O bserved turbulent Mach numbers in e.g., pre sent-day 



proto-plane tary disks are of the order of 0.1-0.5 (Hughes 
et al. 2011). It is typically believed that this disk tur 



bulenc e is driven by the magneto-rotational instability 
(MRI, Balbus fc Hawley 1991). However, in particular 
in the early phases of star formation, disk turbulence 
is possibly driven by the gravitational infall and accre 
tion of gas from the envelope onto the disk 
Hennebelle||2010[ ). The interaction of the MR 



Klessen & 
with self- 
gravitational instabilities may effectively decrea se the ac- 
cretion rate and change the disk morpho logy ( [Fromang 



et al. 2004), also in the primordial case ( [Silk fc Langer 
^Un5| ^. A spectacular example of dynamically important 



magnetic fields in the accretion disks of young stars is the 
generation of high-speed, bipolar jets, which are launched 
by a tower of magnetic pressure and/or b y a centrifugal 
magn etic field (e .g., Blandford & Payne 1982; Pudritz^ 
fc NoSIE I|[T983l ^ontopoulos,.1995, Lynden-BeU 200:^ 



Machida et al.||2006[ [Banerjee fc Pudritz] |2006[ |2007 
Dutfin fc_Pudritz||2009[ for theoretical work on jets anc 



outflows 1^ However, even before that stage, the tur- 
bulent dynamo may have amplified the magnetic field 
enough that it could reduce or even suppress fragmenta- 

^ See also |Beuther et al.| l |2010| | for a recent observational study 
on IRAS 18USy-iY3^, m which -as in our simulations- the turbu- 
lent energies are still dominating the magnetic energy inside the 
dense core. 
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2008 Hennebelle & Teyssier 


2008 


Biirzle et al.) 2010 


Peters et al. 


201ip, thus potentially 



mordial and contemporary star formation. 

We thank Sebastien Fromang for helpful discussions on 
the spectral analysis of non-periodic datasets. C.F. has 
received funding from the European Research Coun- 
cil under the European Community's Seventh Frame- 
work Programme (FP7/2007-2013 Grant Agreement 
no. 247060) for the research presented in this work. 
C.F., R.B., and R.S.K. acknowledge subsidies from the 
Baden- Wiirttemberg-Stiftung (grant P-LS-SPII/18) and 
from the German Bundesministerium fiir Bildung und 
Forschung via the ASTRONET project STAR FOR- 
MAT (grant 05A09VHA). S.S. thanks the German Sci- 



ence Foundation (DFG) for financial support via the 
priority program 1177 "Witnesses of Cosmic History: 
Formation and Evolution of Black Holes, Galaxies and 
their Environment" (grant KL 1358/10). D.S. thanks 
for funding from the European Community's Seventh 
Framework Programme (FP7/2007-2013) under grant 
agreement No 229517. R.B. is funded by the Emmy- 
Noether grant (DFG) BA 3706. Supercomputing time at 
the Leibniz Rechenzentrum (projects pr321o and hl221) 
and the Forschungszentrum Jiilich (projects hhd20 and 
hhdl4) are gratefully acknowledged. The software used 
in this work was in part developed by the DOE-supported 
ASC / Alliance Center for Astrophysical Thermonuclear 
Flashes at the University of Chicago. Figure [3] was pro- 
duced with the open-source visualization software visit. 



REFERENCES 



Abel, T., Bryan, G. L., & Norman, M. L. 2002, Science, 295, 93 
Arshakian, T. G., Beck, R., Krause, M., & Sokoloff, D. 2009, 

A&A, 494, 21 
Balbus, S. A., & Hawley, J. F. 1991, ApJ, 376, 214 
Balsara, D. S., Kim, J., Mac Low, M., & Matliews, G. J. 2004, 

ApJ, 617, 339 
Banerjee, R., & Pudritz, R. E. 2006, ApJ, 641, 949 
— . 2007, ApJ, 660, 479 

Banerjee, R. , Vazquez-Semadeni, E., Hennebelle, P., & Klessen, 

R. S. 2009, MNRAS, 398, 1082 
Bate, M. R., Bonnell, I. A., & Price, N. M. 1995, MNRAS, 277, 

362 

Bate, M. R., & Burkert, A. 1997, MNRAS, 288, 1060 
Begelman, M. C, & Shlosman, I. 2009, ApJ, 702, L5 
Bernet, M. L., Miniati, F., Lilly, S. J., Kronberg, P. P., & 

Dessauges-Zavadsky, M. 2008, Nature, 454, 302 
Beuther, H., Vlemmings, W. H. T., Rao, R., & van der Tak, 

F. F. S. 2010, ApJ, 724, L113 
Biermann, L. 1950, Zeitschrift Naturforschung Teil A, 5, 65 
Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 
Bonnor, W. B. 1956, MNRAS, 116, 351 

Bouchut, F., Klingenberg, C, & Waagan, K. 2007, Numerische 

Mathematik, 108, 7 
— . 2010, Numerische Mathematik, 115, 647 
Brandenburg, A., & Subramanian, K. 2005, Phys. Rep., 417, 1 
Brunt, C. M., Heyer, M. H., & Mac Low, M. 2009, A&A, 504, 883 
Brunt, C. M., Heyer, M. H., Vazquez-Semadeni, E., & Pichardo, 

B. 2003, ApJ, 595, 824 
Burgers, J. M. 1948, Adv. Appl. Mech., 1, 171 
Biirzle, F., Clark, P. C, Stasyszyn, F., Greif, T., Dolag, K., 

Klessen. R. S.. & Nielaba, P. 2010, MNRAS, accepted 

(arXiv:1008.3790) 
Carroll, J. J., frank. A., & Blackman, E. G. 2010, ApJ, 722, 145 
Cho, J., & Vishniac, E. T. 2000, ApJ, 539, 273 
Cho, J., Vishniac, E. T., Beresnyak, A., Lazarian, A., & Ryu, D. 

2009, ApJ, 693, 1449 
Cho, W., & Kim, J. 2011, MNRAS, 410, L8 

Clark, P. C, Glover, S. C. O., & Klessen, R. S. 2008, ApJ, 672, 
757 

Clark, P. C, Glover, S. C. O., Klessen, R. S., & Bromm, V. 2011, 

ApJ, 727, 110 
Contopoulos, J. 1995, ApJ, 450, 616 

Dolag, K., Kachelriess, M., Ostapchenko, S., & Tomas, R. 2011, 
ApJ, 727, L4 

Duffin, D. F., & Pudritz, R. E. 2009, ApJ, 706, L46 
Ebert, R. 1955, Zeitschrift fiir Astrophysik, 37, 217 
Elmegreen, B. G., & Burkort, A. 2010, ApJ, 712, 294 
Elmegreen, B. G., & Scalo, J. 2004, ARA&A, 42, 211 
Federrath, C, Banerjee, R., Clark, P. C, & Klessen, R. S. 2010a, 
ApJ, 713, 269 

Federrath, C, Glover, S. C. O., Klessen, R. S., & Schmidt, W. 

2008a, Phys. Scr. T, 132, 014025 
Federrath, C, Klessen, R. S., & Schmidt, W. 2008b, ApJ, 688, 

L79 

Federrath, C, Roman-Duval, J., Klessen, R. S., Schmidt, W., & 

Mac Low, M. 2010b, A&A, 512, A81 
Field, G. B., Blackman, E. G., & Keto, E. R. 2008, MNRAS, 385, 

181 

Fleck, Jr., R. C. 1983, ApJ, 272, L45 



Folini, D., Walder, R., & Favre, J. M. 2010, in Astronomical 
Society of the Pacific Conference Series, Vol. 429, Astronomical 
Society of the Pacific Conference Series, ed. N. V. Pogorelov, 
E. Audit, & G. P. Zank, 9 

Frisch, U. 1995, Turbulence, the legacy of A. N. Kolmogorov 
(Cambridge Univ. Press) 

Fromang, S., Balbus, S. A., Terquem, C, & De Villiers, J. 2004, 
ApJ, 616, 364 

Fryxell, B., et al. 2000, ApJS, 131, 273 

Gawryszczak, A. J., Mayer, L., Boley, A., & Tasker, E. 2010, in 

EAS Publications Series, Vol. 42, EAS Publications Series, ed. 

K. Gozdziewski, A. Niedzielski, & J. Schneider, 267-274 
Glover, S. C. O., & Savin, D. W. 2009, MNRAS, 393, 911 
Greif, T. H., Johnson, J. L., Klessen, R. S., & Bromm. V. 2008, 

MNRAS, 387, 1021 
Haugen, N. E., Brandenburg, A., & Dobler, W. 2004, 

Phys. Rev. E, 70, 016308 
Heitsch, F., Mac Low, M.-M., & Klessen, R. S. 2001, ApJ, 547, 

280 

Heitsch, F., Slyz, A. D., Devriendt, J. E. G., Hartmann, L. W., & 

Burkert, A. 2006, ApJ, 648, 1052 
Hennebelle, P., Banerjee, R., Vazquez-Semadeni, E., Klessen, 

R. S., & Audit, E. 2008, A&A, 486, L43 
Hennebelle, P., & Chabrier, G. 2008, ApJ, 684, 395 
— . 2009, ApJ, 702, 1428 

Hennebelle, P., & Teyssier, R. 2008, A&A, 477, 25 
Heyer, M. H., & Brunt, C. M. 2004, ApJ, 615, L45 
Hoyle, F. 1953, ApJ, 118, 513 

Hughes, A. M., Wilner, D. J., Andrews, S. M., Qi, C, & 

Hogerheijde, M. R. 2011, ApJ, 727, 85 
Kainulainen, J., Beuther, H., Kenning, T., & Plume, R. 2009, 

A&A, 508, L35 
Kazantsev, A. P. 1968, Soviet Journal of Experimental and 

Theoretical Physics, 26, 1031 
Kitsionas, S., et al. 2009, A&A, 508, 541 
Klessen, R. S. 2000, ApJ, 535, 869 
— . 2001, ApJ, 556, 837 

Klessen, R. S., Heitsch, F., & Mac Low, M.-M. 2000, ApJ, 535, 
887 

Klessen, R. S., & Hennebelle, P. 2010, A&A, 520, A17 
Klingenberg, C., Schmidt, W., & Waagan, K. 2007, Journal of 

Computational Physics, 227, 12 
Kolmogorov, A. N. 1941, Dokl. Akad. Nauk SSSR, 32, 16 
Kritsuk, A. G., Norman, M. L., Padoan, P., & Wagner, R. 2007, 

ApJ, 665, 416 

Kritsuk, A. G., Norman, M. L., & Wagner, R. 2011, ApJ. 727, L20 
Krumholz, M. R., & McKee, C. F. 2005, ApJ, 630, 250 
Krumholz, M. R., McKee, C. F., & Klein, R. I. 2004, ApJ, 611, 
399 

Krumholz, M. R., McKee, C. F., & Tumlinson, J. 2009, ApJ, 699, 
850 

Kulsrud, R. M., Gen, R., Ostriker, J. P., & Ryu, D. 1997, ApJ, 
480, 481 

Larson, R. B. 1969, MNRAS, 145, 271 
— . 1981, MNRAS, 194, 809 

Lazar, M., Schlickeiser, R., Wielebinski, R., & Poedts, S. 2009, 

ApJ, 693, 1133 
Lemaster, M. N., & Stone, J. M. 2008, ApJ, 682, L97 
Li, Y., Klessen, R. S., & Mac Low, M.-M. 2003, ApJ, 592, 975 
Li, Z., Wang, P., Abel, T., & Nakamura, F. 2010, ApJ, 720, L26 
Lynden-Bell, D. 2003, MNRAS, 341, 1360 



16 



Federrath, Sur, Schleicher, Banerjee, Klessen 



Mac Low, M.-M., & Klessen, R. S. 2004, Reviews of Modern 
Physics, 76, 125 

Machida, M. N., Matsumoto, T., & Inutsuka, S. 2008, ApJ, 685, 
690 

Machida, M. N., Omukai, K., Matsumoto, T., & Inutsuka, S. 

2006, ApJ, 647, LI 
Maki, H., & Susa, H. 2004, ApJ, 609, 467 
McKee, C. F., & Ostriker, E. C. 2007, ARA&A, 45, 565 
Miniati, F., & Bell, A. R. 2011, ApJ, 729, 73 
Nakamura, F., & Li, Z. 2008, ApJ, 687, 354 
Odgers, G. J., & Stewart, R. W. 1958, Reviews of Modern 

Physics, 30, 1017 
Omukai, K., Tsuribe, T., Schneider, R., & Ferrara, A. 2005, ApJ, 

626, 627 

O'Shea, B. W., & Norman, M. L. 2007, ApJ, 654, 66 
Ossenkopf, V., Krips, M., & Stutzki, J. 2008a, A&A, 485, 917 
— . 2008b, A&A, 485, 719 

Ossenkopf, V., & Mac Low, M.-M. 2002, A&A, 390, 307 
Padoan, P., & Nordlund, A. 2002, ApJ, 576, 870 
Padoan, P., Nordlund, A., & Jones, B. J. T. 1997, MNRAS, 288, 
145 

Passot, T., & Vazqucz-Semadeni, E. 1998, Phys. Rev. E, 58, 4501 
Pension, M. V. 1969. MNRAS, 144, 425 

Peters, T., Banerjee, R., Klessen, R. S., & Mac Low, M. 2011, 
ApJ, 729, 72 

Price, D. J., & Bate, M. R. 2007, MNRAS, 377, 77 

Price, D. J., & Federrath, C. 2010, MNRAS, 406, 1659 

Price, D. J., Federrath, C, & Brunt, C. M. 2011, ApJ, 727, L21 

Pudritz, R. E., & Norman, C. A. 1983, ApJ, 274, 677 

Pudritz, R. E., & Silk, J. 1989, ApJ, 342, 650 

Ryu, D., Kang, H., Cho, J., & Das, S. 2008, Science, 320, 909 

Scalo, J. M., &; Pumphrey, W. A. 1982, ApJ, 258, L29 

Schekochihin, A. A., Cowley, S. C, Taylor, S. F., Maron, J. L., & 

McWilliams, J. C. 2004, ApJ, 612, 276 
Schleicher, D. R. G., Banerjee, R., Sur, S., Arshakian, T. G., 

Klessen, R. S., Beck, R., & Spaans, M. 2010, A&A, 522, A115 



Schleicher, D. R. G., GaUi, D., Glover, S. C. O., Banerjee, R., 
Palla, F., Schneider, R., & Klessen, R. S. 2009, ApJ, 703, 1096 

Schmidt, W., Federrath, C., Hupp, M., Kern, S., & Niemeyer, 
J. C. 2009, A&A, 494, 127 

Schneider, N., et al. 2010, A&A, accepted l|arXiv;1001.2453|l 

Silk, J., & Langer, M. 2006, MNRAS, 371, 444- 

Solomon. P. M., Rivolo, A. R., Barrett, J., & Yahil, A. 1987, ApJ, 
319, 730 

Stutzki, J., Bensch, F., Heithausen, A., Ossenkopf, V., & 

Zielinsky, M. 1998, A&A, 336, 697 
Sur, S., Schleicher, D. R. G., Banerjee, R., Federrath, C., & 

Klessen, R. S. 2010, ApJ, 721, L134 
Tan, J. C., & Blackman, E. G. 2004, ApJ, 603, 401 
Tassis, K. 2007, MNRAS, 382, 1317 

Tavecchio, F., Ghisellini. G., Foschini, L., Bonnoli, G., Ghirlanda, 

G., & Coppi, P. 2010, MNRAS, 406, L70 
Truelove, J. K., Klein, R. I., McKee, C. F., HoUiman, II, J. H., 

Howell, L. H., & Greenough, J. A. 1997, ApJ, 489, L179 

Vazquez-Semadeni, E. 1994, ApJ, 423, 681 

Vazquez-Semadeni, E. 2010, lAUS, in press (arXiv:l6l0.5157) 
Vazquez-Semadeni, E., Ryu, D., Passot, T., Gonzalez, il. t ., & 

Gazol, A. 2006, ApJ, 643, 245 
Vishniac, E. T. 1994, ApJ, 428, 186 

Waagan, K. 2009, Journal of Computational Physics, 228, 8609 
Waagan, K., Federrath, C., & Klingenberg, C. 2011. Journal of 

Computational Physics, in press ( arXiv:1101.3007 ) 
Wada, K., Meurer, G., & Norman, C. A. 2002, ApJ, 577, 197 
Walder, R., & FoUni, D. 1996, A&A, 315, 265 
Wise, J. H., & Abel, T. 2007, ApJ, 665, 899 

Xu, H., Li, H., Collins. D. C., Li, S., & Norman, M. L. 2009, ApJ, 

698, L14 
— . 2010, ApJ, 725, 2152 

Yoshida, N., Omukai, K., & Hernquist, L. 2008, Science, 321, 669 



APPENDIX 

POWER-LAW DENSITY PROBABILITY DISTRIBUTION FUNCTIONS 
Given a power-law, spherical density distribution with exponent a, 



P{r) = Po 



ro 



(Al) 



we can derive the power-law exponent of the probability distribution function (PDF) of the density. The differential 
volume of the spherical gas distribution is a shell with 



dV = Anr dr . 



The volume-weighted PDF of density, Pp is given by 



dV dV dr 
dp dr dp \po 



-(^1+1) 



(A2) 

(A3) 
Li 



The PDF, of the logarithmic d ensity, s = \n(p/{p)), is related to the PDF, pp, of the linear density (see, e.g 
et al.||2003[ [Federrath etar]|2008b[ ) by 

Ps^ PPp , (A4) 

such that 



Ps oc 



cx exp 



3s 



(A5) 



REYNOLDS NUMBER DEPENDENCE OF THE DYNAMO GROWTH RATE 

In the regime of intermediate to large magnetic Prandtl number, Pm = v/rj > 1, where v is the kinematic viscosity 
and r] is the magnetic diffusivity, the dynamo growth rate, is determined by the vorticity, a; = |V x u| on small 
scales. We thus make the ansatz 

n = au}, (Bl) 

with a numerical constant a of order unity, which is determined by fitting the growth rate in a time evolution plot of 
the magnetic field, from numerical experiments (see, e.g., Fig. From lyu^ — fe^u^/^int with another dimen sion less 



constant b, the integral scale £i^t of the turbulence and the velocity dispersion ct„ on that scale, and equation (Bl I, it 
follows that the growth rate 

(B2) 
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Fig. 13. — Test on the numerical diffusion of the magnetic field in the HLL3R scheme for ideal MHD ( [Waagan et al.||2011[ ). Three 
runs with Jeans resolutions of 8, 16, and 32 cells are shown. For clarity, the curves for 16 and 8 cells were smtted down by iU and 20%, 
respectively. The initial decrease of B/p^^^ is due to the anisotropic geometry of the initial field structure. At late times during the 
collapse, tfie field becomes approximately isotropic due to the bending of the field lines to an hour-glass shape. In this regime, we can 
estimate upper limits for the effects of numerical diffusion, indicated by the dashed lines for p > 10~^ gcm~^. 



with the eddy turnover time ^od = ^int/2cri,. This shows that the growth rate for dynamo amplification is expected 
to increase with the square root of the Reynolds number Re = cr^^int/i^- For magnetic Prandtl numbe r, Pm = 1, the 
kinematic and magnetic Reynolds numbers are the same: Re = Rm, such that fl oc Rm^^^ (see alsq^ Haugen et al. 
2004 1 . Since both i/ and r] are controlled by the numerical cutoff scale in our simulations (see section |2.2p, Pm 1. 



TEST ON THE NUMERICAL DIFFUSION OF THE HLL3R IDEAL MHD SCHEME 



Here we test the numerical diffusion of the magnetic field in the new HLL3R scheme for ideal MHD ( Waagan 
et al. [2011), implemented in the FLASH code. We use the same setup as in section [2] but set all initial velocities to 
zero. Unlike the turbulent magnetic field structure in the production runs, we start with a uniform magnetic field in 
z-direction with B = 10~^ G. The magnetic field is dragged in and develops an hour-glass-shaped magnetic field. 
In Figure 



13 



we show the magnetic field strength divided by p^^^ as a function of the density, measured within 
the collapsing, central Jeans volume (see analysis described in section 2.3 1. The magnetic field grows with increasing 
density, but initially, it grows shallower than B cx p^/^ due to the geometry of the field. Only for a fully isotropic field, 
B (X p^/^ for ideal MHD. Here, however, the magnetic field is initially very anisotropic (containing a z-component 
only), and thus, the magnetic field is initially only compressed along the x—, and y-directions, while the density 
increases along all three axes. Due to the bending of the field lines, the field becomes more isotropic in the late stages 
of the collapse. Figure 13 indeed shows a clear flattening of Bj p^l^ at high densities, when due to the bending, the 



magnetic field has developed roughly similar magnetic field strengths in all three dir ections. Thus, only then, we expect 
B (X p^l^ for ideal MHD. Power-law fits in this regime are indicated in Figure |l3| The curves are still falling with a 
roughly constant dependence on density, giving an upper limit on the effects of numerical diffusion. Three runs with 
Jeans resolutions of 8, 16, 32 cells were analyzed in this diffusion test. We find a departure from perfect flux-freezing, 
B cx with e k, 0.005, 0.003, and 0.002 for a Jeans resolution of 8, 16, and 32 cells, respectively. Thus, the effects 

of numerical diffusion are resolution-dependent, as expected, however, the level of numerical diffusion of the magnetic 
field with respect to ideal MHD is generally quite small. 

In this simple test, we only studied a particular situation, motivated by the density and magnetic field regimes 
analyzed in the main part of the paper. It should be noted, however, that the properties of numerical diffusion depend 
on the physical and numerical parameters of the modeled system. 



